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Abstract 

Bicontinuous, interfacially jammed emulsion gels (bijels) represent a class of soft 
solid materials in which interpenetrating domains of two immiscible fluids are stabi- 
lized by an interfacial colloidal monolayer. Such structures can be formed by arrested 
spinodal decomposition from an initially single-phase colloidal suspension. Here we 
explore by lattice Boltzmann simulation the possible effects of using magnetic colloids 
in bijels. This may allow additional control over the structure, during or after forma- 
tion, by application of a magnetic field or field gradient. These effects are modest for 
typical parameters based on the magnetic nanoparticles used in conventional ferroflu- 
ids, although significantly larger particles might be appropriate here. Field gradient 
effects, which are cumulative across a sample, could then allow a route for controlled 
breakdown of bijels as they do for particle-stabilized droplet emulsions. 

1 Introduction 

The existence of particle-stabilized emulsions has been known for over a century [I] . Recently, 
understanding of their formation has advanced considerably [2], and stimulated a wider 
interest in using liquid-liquid interfaces to direct the self-assembly of nanoparticles (3j HJ [5] . 
Recent simulation work in our group [6j [7] has addressed a new route to directed self assembly 
of colloids (including nanocolloids), creating ffuid-bicontinuous structures. If the interfacial 
particles solidify by jamming in such a case, it should impart macroscopic rigidity in three 
dimensions [HI IS]. The resulting 'bijel' (bicontinuous interfacially jammed emulsion gel), 
comprising a solid matrix permeated by a pair of bicontinuous fluids, could have potential 
applications, for instance as a 'membrane contactor' for catalytic applications [6| |9| [TO]. 

The basic protocol for making bijels starts with a sample in the single phase region of a 
binary fluid mixture, in which are suspended colloidal particles having roughly equal affinity 
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for the two fluids. There follows a quench in which the fluids demix, sweeping up particles 
to the interface, where they adopt a contact angle close to 90° (neutral wetting). As the 
interface reduces its area by the usual coarsening process, the colloids become jammed into 
close proximity. There is then an energy barrier cue, where e = ana 2 , to detachment of a 
particle of radius a from a fluid-fluid interface of tension a, where a is a geometrical parameter 
dependent on the local particle packing. (Naively one might expect a ~ 1, although we have 
presented numerical evidence for smaller values [7]). But unless a is very small indeed, the 
detachment barrier remains large compared to /c#T for colloids larger than a few nm [2], in 
which case the jammed bijel structure should be permanently arrested. Such fully arrested 
bijels have recently been created experimentally [HI [El EES], and do exhibit a yield stress. 
In practice the colloidal particles used - although repulsive at moderate separations - might 
be pushed by the interfacial forces into a primary flocculation minimum creating quasi- 
permanent bonds [13]. However, this takes several hours (at least in the systems studied in 
[T5] ) and so can be viewed as a final lock-in of a structure whose formation is governed by 
the interplay of interfacial and repulsive colloidal forces. 

In this paper we present new simulations addressing the case where the colloidal parti- 
cles have magnetic (or dipolar) interactions as well as the short-range repulsion considered 
previously. Magnetic colloids are of course well known experimentally as the basis of fer- 
rofluids - liquids whose material properties, particularly viscosity, are strongly sensitive to 
external magnetic fields [2]. A very strong field-dependence has also been demonstrated 
for emulsion droplets stabilized by magnetic particles [15] : specifically, it was found that 
for large (millimetre) emulsion droplets, the field gradient of an ordinary bar magnet was 
enough to pull magnetic particles off the fluid-fluid interface, breaking the emulsion. That 
work involved superparamagnetic particles rather than the permanent dipoles considered in 
our simulations, but the fields were high enough to give near-saturation of the dipole mo- 
ments [15] , and in this limit the two systems should behave similarly. This suggests a wider, 
and so far unexploited, avenue for using magnetic particles at fluid-fluid interfaces to create 
field- sensitive, self-assembled materials. We explore this here in the bijel context. 

The rest of this paper is organized as follows. In Section [2J we give a brief account of the 
simulation methodology; this combines elements from our earlier work on bijels [7J and on 
magnetic colloids |16] . We then quantify in Section [3] the interfacial morphologies of bijels 
and droplets containing magnetic particles, created in the presence or absence of a uniform 
magnetic field. In Section H] we consider cases where a field or field gradient is applied late in 
the formation of a bijel, and investigate the extent of the resulting structural perturbations. 
Section \5\ contains a discussion. 

2 Simulation methods 

We use the lattice Boltzmann (LB) method [TTJ for a binary fluid [TS] incorporating spherical 
solid particles [2D], modifying a standard 'bounce-back on links' method fTDJ [2JJ, to allow 
for the presence of a binary solvent [20, 22J. The two solid-fluid interfacial tensions are 
exactly equal, with the interfacial thermodynamics implemented as reported previously [20] . 
We choose for simplicity a pair of fluids with equal density p and viscosity t]. Their phase 
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diagram is also symmetric, and described by the free energy functional [23] 



F[<0] = dV (Ai) 2 j2 + B^j\ + k( W) V 2 ) 



(1) 



where the order parameter ip describes the fluid composition, and the model parameters 
A < 0, B, and k control the fluid-fluid interfacial tension a and thickness £ [23J. The ip = 
isosurface defines the position of the fluid-fluid interface. 

The binary fluid is initialized to be well mixed and at rest, with a small amplitude random 
noise added to the ip field to induce spinodal decomposition. The colloids are initially 
positioned at rest randomly throughout the system: in effect therefore we quench from 
infinite temperature. We choose a deep quench in the sense that thermal fluctuations of i/j 
remain negligible: without particles, the phase separation proceeds purely by amplification of 
the initial noise. Time- dependent, thermal noise is however fully included in the description 
of fluid momentum, using a method reported previously [21]. As a result, the colloids undergo 
realistic, many-body Brownian motion. 

Unless otherwise stated below, our parameters are set as follows. First, we choose —A = 
B = 0.002, and k = 0.0014, giving an interfacial thickness of £ = 1.14 and tension a = 
1.58 x 10 -3 ; we set the (scalable) fluid density p — 1, and choose viscosity r\ — 0.1. Here and 
below, all physical variables are expressed in lattice units [23] , except where SI units are given 
explicitly. For numerical efficiency, the colloid radius a is made as small as is compatible 
with acceptably accurate simulation [BJ [23 [21]; we choose a = 2.3, so that e = 0.026. Our 
default noise setting is k^T = 2.13 x 10~ 5 , so that e/ksT = 1230. We can choose to interpret 
these parameters as representing a short-chain hydrocarbon/water or hydrocarbon/alcohol 
mixture, with laboratory parameters p = 1000 kg m" 3 , 77 = 9.3 x 10~ 4 N m _2 s (close 
to that of water at ~ 300K) and a = 6.1 x 10~ 2 N m _1 . This corresponds to a physical 
particle radius of a = 5.1 nm. We have in addition performed several runs at a temperature 
fcgT = 2 x 10~ 4 , approximately ten times larger than the default value. This allows us 
effectively to probe longer timescales (since the Brownian diffusivity increases by the same 
factor) while maintaining a value of the capillary energy that is still very large compared to 
thermal energies [ejk^T ~ 130). 

The volume fraction of colloids is <fi = 0.20 in all simulations reported here. However, 
in the bijel these particles get packed into close proximity forming a network of repulsive 
contacts. LB does not properly resolve the lubrication films at these contacts; to correct 
for this we introduce a normal lubrication force for surface separations h = r — 2a < h c = 
0.3. Such lubrication forces slow the code down considerably once there are large clusters 
of particles in mutual lubrication contact; to minimize this, a soft-core thermodynamic 
repulsion is introduced so that interfacial particles maintain a separation somewhat larger 
than 2a [6]. (Such a repulsion is often present physically, e.g., for charge-stabilized colloids.) 
The chosen 'soft-core' short-range potential U sc (h) as a function of the particle surface 
separation h = r — 2a obeys 



where u(h) = , y(ho/h)' / , and h c is a cut-off separation. (Thus U sc is a 'cut-and-shift' power 
law potential which vanishes smoothly at h c .) The parameter set is chosen as 7 = lO/c^T, 
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ho = 0.1, v = 1.0 and h c = 0.25, as benchmarked previously [7J [16]. All simulations have 
been done on a D3Q19 lattice with the system volume A 3 = 64 3 and with either periodic 
or fixed-wall boundary conditions. For the latter conditions, a normal lubrication force is 
introduced for colloids at distances h < hi u b = 0.5 between the plane of the wall and the 
surface of the colloid; a neutral wetting condition is maintained between such boundary walls 
and the fluid-fluid interfaces. 

A dipole in an external magnetic field B is subject to a force F, = mjSj.VB and a torque 
Tj = rriiSi x B, where rrii is the magnitude of dipole moment i and Sj its direction. (We take 
rrii = m for all particles.) The dipole-dipole interaction is given directly as 
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where the prefactor fio/4-K = 10 7 H/m. The relative importance of the dipole-dipole inter- 
action, compared to thermal energies, is characterized by the dimensionless parameter 

= Ho rn 2 

An8a 3 k B T' 1 ; 

For ferrofluid applications, A values in the range 2 — 8 are readily available but larger values 
remain relatively uncommon (see [16] and references therein). To calculate the long-range 
dipolar interactions for a periodic simulation box, Ewald summation is used following the 
procedure and parameter choices laid out in |TB]: this is unaffected by the switch from a 
single solvent in that paper to a binary one here. In the case of fixed wall boundary conditions 
the long-range dipolar interaction was directly calculated from Eqj3j 

Below we examine the behavior of bijels stabilised by magnetic colloids under either a 
uniform magnetic field or a field gradient. A uniform field B directly creates torque but 
not net forces on the dipoles; by aligning the dipoles, it alters the dipole-dipole interactions. 
A uniform field applied within the simulation box clearly obeys V.Bo = 0, and can be used 
with either periodic or fixed- wall boundary conditions. In the case of a field gradient, we 
restrict attention to fixed- wall boundary conditions since the requirement of V.Bo = is 
then not easily satisfied in the periodic case. The field pattern chosen is in fact that of a 
point magnetic dipole (much larger in magnitude than those on the particles) held at a fixed 
orientation at a location outside the simulation box. 

Most of our simulations have been evolved for 5 x 10 5 timesteps; this is on the microsecond 
timescale for nanocolloids [6] and is long enough to reach a near-steady state with our 
chosen lattice size (A 3 = 64 3 ), at least so long as one disregards the slow aging effects 
exploreded in [7]. For one such run, a single Intel 2.4GHz processor requires 270 hours. 
Parallel computation using 8 cores of a cluster of 3GHz Intel Dual-core processors reduce 
that computational time to about 26 hours. 



3 Fluid and magnetic morphologies 

In this Section we present the morphologies of fluid domains arising under various condi- 
tions. We vary the initial volume ratio of fluids (ipo), the dipolar coupling constant (A), and 
the strength of a uniform external field (B ). We always use the same initial conditions 
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with a uniformly mixed fluid phase and random positions of colloids. We use the default 
temperature value except where otherwise stated. 

3.1 Bicontinuous phase and droplet phase 

The fluid domain morphologies resulting from our chosen quench protocol have been pre- 
sented elsewhere for nonmagnetic colloids (A = 0) [7j. By varying the composition variable 
ipo away from the symmetric value ipo = 0, the percolation threshold ip p between bicontinuous 
and droplet morphologies was estimated to lie within the range from 0.3 to 0.4. 

Figure [I] presents snapshots for the domain structure at the end of the simulation with 
a dipolar coupling constant A = 4 for both a symmetric and a strongly asymmetric (ipo = 
0.4 > ip p ) quench, corresponding to fluid ratios of 50:50 and 70:30 respectively. (We use 
periodic boundary conditions; results for closed-box boundary conditions are similar.) In 
each case the basic morphology is broadly similar to that found for A = in previous 
work |7|. This is not surprising since, although the magnetic interaction is several times 
larger than &bT, the balance between capillary forces and hard-core repulsions involves a 
much higher energy scale than this (e = 1230k bT). Within the particle layer, however, the 
magnetic interactions promote the dipolar particles to align in a nose-to-tail fashion in all 
cases; this requires the dipole vectors to lie tangent to the local monolayer. The curvature of 
this monolayer then prevents the formation of a uniform monodomain of particles, creating 
frustration. Thus in the absence of an external field we see no long range order and expect 
zero mean magnetization (although in practice a small nonzero value is measured due to 
finite size effects). In FigJTJ we identify a representative droplet within the depercolated 
regime; in Section I3~21 we study more carefully the local magnetic ordering at the surface of 
this droplet. 

We have tried increasing A up to A = 40, well above the usual experimental range for 
magnetic colloids. We find that for symmetric quenches, although the magnetic ordering 
within the particle layer increases (longer chains of head-to-tail contacts) the bicontinuous 
morphology of fluid domains is not significantly altered. The same applies for the asymmetric 
quench at ipo = 0.3 whose structure can be seen in FigJ^a) below, where it represents an 
initial state prior to switching on a magnetic field. 

3.2 Droplet covered by magnetic particles 

Figure |5] shows a droplet chosen from the snapshot in Figfjjb), where the magnetic particles 
on its surface have aligned to reduce their interaction energies. As discussed in ref. [7], 
this droplet is aspherical. Figj2]^a) shows the droplet with different viewpoints rotating by 
a full turn counter-clockwise; hence pictures 1 and 6 are the same. The particles achieve a 
local packing that is mostly close to hexagonal ordering, with dipolar orientations primarily 
nose-to-tail, following the curvature of the surface. In pictures 1 and 3, domains of parallel 
magnetic alignment are seen, but with opposite orientations on the two faces of the droplet. 
The intermediate viewpoints (2 and 4) show nontrivial defect structures, as do Figj2]^b,c) 
showing the top and bottom of the droplet. This arrangement of dipoles shows considerably 
more frustation than strictly necessary to accommodate a locally aligned vector field on a 
surface of spherical topology (see, e.g. ref. [2Z])- Although the observed magnetic alignment 
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(a) 




(b) 

Figure 1: The morphologies in a system with periodic boundary conditions for A = 4 at 
k B T = 2.133 x 10~ 5 ; system size is A 3 = 64 3 . The right panels shows the fluid- fluid interfaces 
plus colloids represented as translucent spheres of radius a. Within each such sphere, a rod 
is shown to represent the magnetic dipole whose poles are painted red and white. The left 
panels show the same images without the particles to enable fluid domains to be seen more 
clearly, (a) Symmetric quench: ipo = 0.0. (b) Asymmetric quench: ipo = 0.4. One droplet 
(circled in red) is chosen to see the geometry of magnetic particles on the surface, discussed 
in Section l3~2l 
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might anneal somewhat if the simulation were run for much longer, at least some of these 
defect structures are likely to be permanently trapped as a result of the underlying defects 
in the particle packing. Note that for a perfect hexagonal lattice, there is a magnetic energy 
gain of 4XksT per particle in creating a perfect in-plane head-to-tail alignment compared 
to one in which the dipoles point at 30 degrees (but still in plane) to the lattice vector; but 
with a defect-ridden packing (locked in by capillary forces) local as well as global frustration 
may prevent this energy gain from being realized. 




(c) (d) (e) 

Figure 2: Dipolar ordering on the surface of the droplet, (a) Snapshots for sides of the 
droplet rotating counter-clockwise, (b) Viewed from above, (c) Viewed from below, (d) 
The directions of dipoles covering the droplet are shown with particles coloured according 
to Uf/ksT. A frustrated particle (orange) is found among dipoles on the surface, having 
dipolar energy U d /k B T = -4.70. (e) The local dipolar energy distribution of surface colloids. 

In general, the dipolar energy for an individual particle can be defined from EqJ3] as 
Uf = J2j Ufj, where to maintain consistency with periodic boundary conditions, the sum is 



7 



restricted to < A/2. In the Ewald method the remaining contribution is computed in 
Fourier space, but in the absence of long range order this term should not dominate, and we 
can usefully study the distribution of local energies with it neglected. We show in Figf2]^d) 
an image of the droplet with particles color-coded according to their magnetic energy, and 
in FigfS^e) a histogram of the dipolar energies observed. The single most frustrated particle 
(full red, bottom right in (d)) is in fact not attached (or barely so) to the fluid-fluid interface. 
(The method used to identify particles includes all those whose centres lie within a certain 
distance of the interface, but depending on the arrangement of other particles this may 
or may not signify attachment.) Of the fully attached particles the most frustrated is at 
upper left centre (translucent orange) with an energy of — 4.7/cbT. This is surrounded by a 
hexagonal ring on which the dipolar axis rotates through a full turn. 

3.3 Effect of uniform external field 

In the absence of an applied external field, magnetic dipoles have only local orientational 
correlations. However, under the external field, dipole moments preferentially align along 
the field direction. At very high field strength, completely aligned dipoles can be expected. 
At intermediate values, the magnetic saturation for noninteracting permanent dipoles is 
predicted by Langevin theory [14] as a function of = mB/hsT. A recent study of 
ferrofluids in a uniform field [26] has shown that the noninteracting Langevin theory is still 
applicable in systems with moderate dipole-dipole interactions. 

Figure |3] shows domain morphology snapshots for simulations in which the quench was 
made with a uniform magnetic field present from the outset, with as = 2 and 20 respectively 
and with A = 4. To speed approach to the final, locally equilibrated structure, we performed 
these runs with enhanced diffusion (achieved by setting fcgT = 2 x 10~ 4 ). The applied uni- 
form field is given as B = (0, 0, B z ) and the corresponding Langevin values for the magnetic 
saturation J2i{si)z/N are s = 0.54 and s = 0.95. The measured values are instead s = 0.69 
and s = 0.93, presumably reflecting the combined effects of surface confinement and dipolar 
interactions. At as = 2, the initial domain morphology is only weakly perturbed by the 
interactions of dipoles either with each other, or with the external field. At the higher field 
strength, as = 20 corresponding to B = 0.31T in the laboratory, significant distortion of 
the fluid-fluid interface is seen, stretching it along the field direction. Although the capillary 
and hard-core forces remain strong (we now have e/ksT ~ 130 for the capillary energy), 
at these field levels the requirement of dipolar alignment with the field, combined with the 
energetic preference for head-to-tail dipolar chains, can compete with them. However, the 
degree of distortion remains relatively modest, and, at least for a symmetric quench, fluid 
bicontinuity is maintained. 

3.4 Time evolution of characteristic domain lengths 

The morphology of the bijel structure in the presence of a field can be characterized via length 
scales L Xi y >z defined J2S] as ^L x = l/d xx , etc., where d a/3 = S r ^?^( r )^f VK r )/ 2 r V^ 1 ")- Here 
d D is a symmetric discrete derivative and the sum is over all fluid sites. These characteristic 
lengths measure in essence the mean distance between interfaces along a line drawn along 
each of the coordinate axes; with a field along z, directions x, y should be statistically 
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Figure 3: Final morphologies for A = 4 at ksT = 2 x 10 4 . System size is A 3 = 64 3 ; periodic 
boundary conditions are used, (a) In zero field, (b) With dimensionless field strength a B = 2. 
(c) With a B = 20. 
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equivalent. (However, once the domain size becomes comparable to the simulation box, their 
L values may differ significantly in any instance.) For the case of zero field, the structure is 
isotropic and a more convenient measure of the length scale is L(t) = 2tt J S(k)dk/ J kS{k)dk 
where S{k) is the usual density- density correlator [23] . 

Figure Ufa) presents the time evolution of L(t) for symmetric and asymmetric quenches 
in the absence of a field, with A = 4 and, for comparison, A = 0. At the chosen (higher) 
temperature, the correlated relaxation of dipolar particles leads to a decrease in the do- 
main length scale L(t) at given time, compared to the case without dipoles. (At the default 
temperature, not shown, there is instead a modest increase.) Fig. H^b-d) show L xyz for sym- 
metric quenches in uniform fields of strength as = 2, 20 again with A = 4. As visible already 
in Figj3j whereas for as = 2 any elongation lies within the statistical scatter arising from 
a near-isotropic structure, for the higher field we find significant stretching along the field 
direction with a domain length along z about 50% larger than in the transverse directions. 

3.5 Energetic relaxation of magnetic colloids 

Our bijel quenches always start with the colloids positioned at random; in the magnetic 
case this corresponds to a high temperature limit (A = 0) as the initial state. Thus when 
the quench is performed, the fluid phase separation and the aggregation tendency of dipolar 
colloids are switched on concurrently. As in the nonmagnetic case, particles become jammed 
at the fluid-fluid interfaces causing arrest of L(t) and then a slow aging of this quantity. 
With the temperature at its default value, we observe in addition a very slow relaxation 
of the dipole-dipole energy, FigfS^a), suggestive of incomplete equilibration of the magnetic 
degrees of freedom, even under conditions where the particle packing is considered fixed. As 
explained above, the magnetic relaxation can be speeded up by increasing the diffusivity 
using an enhanced temperature (fegT = 2 x 10 4 ) while retaining a large value e ~ 130/cbT 
for the capillary energy scale. Once this is done, the observed behavior is consistent with 
relaxation to a state of local equilibrium for the given particle packing; the dipolar energy 
continues to evolve slowly, but at a rate set by slow aging of the packing itself [7]. Figj5]^b,c) 
show the resulting relaxations of the dipolar energy and also the magnetic field energy per 
particle, in quenches at A = 4 with a# = 2 and 20. The field energy shows a similar aging 
process in a weak field but, at = 20, saturates much sooner. This is not surprising 
as for strong fields the particles will align fully with the field direction regardless of their 
disposition at the interface. The final field energies per particle (with U B = —mJ2i^i-^o) 
are U B /NksT = —1.37 and —18.65 for = 2,20 respectively; these values follow directly 
from the magnetic saturation values s = 0.69 and s = 0.93 quoted previously. 

Table I presents the dipole-dipole and short-range repulsive contributions to the mean 
energy per particle for all runs so far discussed. Both quantities are averaged over a time 
window comprising the last 50, 000 timesteps in each run. This averaging improves the 
statistics, although significant run-to-run variations remain (particularly in the default tem- 
perature samples, which are not fully equilibrated magnetically). This is not surprising since 
the arrested domain length-scales L xyz are of order half the simulation box size, so the system 
is not self-averaging. For magnetic colloids in a single fluid [29], the equilibrium dipolar and 
soft-core energies at A = 4, as = and = 0.20 are measured as U d /ksT = —4.904 ± 0.004 
per particle and U S€ jk^T = 0.0326 ± 0.0001. Values of the latter are higher in magnetic 



10 




le+05 2e+05 3e+05 4e+05 5e+05 




le+05 2e+05 3e+05 4e+05 5e+0f 
t 

(c) 



le+05 2e+05 3e+05 4e+05 5e+05 
t 

(d) 



Figure 4: (a) The characteristic length scale, L(t) for magnetic bijels in zero field. Black 
and red solid lines are ipo = 0.0 and 0.4 respectively, with A = 4. Green and blue dashed 
lines are ipo = 0.0 and 0.4 with A = 0. Frames (b-d) show L X}V)Z , the characteristic length 
scales in an applied uniform field, for dipole-dipole strength A = 4. The black solid line has 
dimensionless field strength ag = 2; the green dashed line has as = 20. For the parameters 
selected (here with /c#T = 2 x 10~ 4 , 40 lattice units corresponds to 88.4nm, while 500,000 
time steps corresponds to 275ns, or about 5 Brownian times (the time taken for the particle 
to diffuse its own radius).) 
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Figure 5: (a) The dipole-dipole energy per particle U d /NksT at A = 4 for ip = 0.0 
(solid line) and ipo = 0.4 (dashed line) with ksT = 2.133 x 10 -5 . This shows incomplete 
relaxation of the magnetic degrees of freedom, as revealed by the remaining two plots for 
which ksT = 2 x 10 -4 , giving enhanced rotational diffusion of the particles without creating 
a significant tendency of colloids to leave the interface [7J. These show (b) U d /iV/cgT and 
(c) U B jNksT at A = 4, if) = 0; the black solid line has as = 2 and the green dashed line 
has as = 20. 

bijels than in conventional ferrofluids, as a result of the jamming of colloids into close con- 
tact at the fluid-fluid interface. Among equilibrated samples, the dipolar energies are made 
somewhat more negative by the application of a field, consistent with a suggestion that the 
field modestly reduces the magnetic frustration of the particles at the interface. 

4 Structure deformation under external fields 

In this Section we will discuss the effect of switching on a uniform field or field gradient on 
a pre-existing magnetic bijel. Recall that a uniform field creates only torques which cause 
dipoles to arrange along the direction of field; a field gradient creates both forces and torques. 
For this reason the two cases are quite different in their effects. 

4.1 Uniform field 

We first examine the effect of switching on a uniform field, on a pre-existing magnetic bijel. 
From an applications viewpoint, it is interesting to see whether a morphological transition, 
for instance from a bipercolating to a droplet or lamellar morphology, can be brought about 
by application of a field. To this end we consider a bicontinuous but asymmetric bijel with 
ipo = 0.3 in a periodic box of A 3 = 64 3 . (This is close to the percolation threshold; one 
might hope to see e.g. a loss of bicontinuity transverse to the applied field.) To maximise 
the chances of seeing such an effect, we choose magnetic colloids with extremely high dipolar 
interactions, A = 40, retaining the usual volume fraction = 0.2. (We have also performed 
runs with lower A; for A < 8 the effect of field on structure is very modest.) To see large 
effects we need also very high field strengths, with a B = mB /kBT in the range 50-1000. 
Although these A and as values are very high, they are not completely unfeasible: for 
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Table 1: Dimensionless energies in zero and uniform fields. 

instance, A = 40 could be achieved using monodomain magnetite particles (magnetization 
density 4.8 x 10 5 Am _1 ) of 14.7nm radius, and for such particles a B = 1000 equates to 
Bq ~ 0.63T. This radius is close to the theoretical size limit for magnetite monodomains [30] , 
beyond which particles become superparamagnetic rather than having permanent dipoles. 
As mentioned previously, however, in the presence of large fields as considered here, larger 
particles will anyway show saturating magnetization and while we have not simulated the 
superparamagnetic case, they should behave similarly to permanent dipoles with very large 
A values so long as the field remains in place. 

Figure [6] shows snapshops of the final structures and graphs of the time evolution of the 
characteristic domain length-scales for a B = 50,100,1000. We have reverted here to the 
default temperature so that e = e/k B T = 1230. Thus the ratio of dipole-dipole to capillary 
energies is A/e = 0.0325, and the ratio of field to capillary energies a B /e lies in the range 
0.04 — 0.81. The field is applied along the z direction which is here horizontal. The length 
scale plots show clear deformation of the structure; this takes some time to establish, but then 
appears to saturate. Even for the largest field we do not see depercolation of the bicontinuous 
structure in the final state, although this might possibly happen in a much larger sample, 
where the visible tendency to form horizontal sheets and/or tubes could eventually close off 
all fluid pathways transverse to the field. 

The structures seen for a B > 100 are quite similar to the one shown in Figj3]^c), for which 
the value of a B /e = 0.15 lies in the same range. This supports the idea that the field- induced 
effects on morphology are governed primarily by this ratio of field energy to capillary energy. 
Qualitatively, it does not seem crucial whether the field is applied during the initial formation 
of the bijel or switched on afterwards. To check whether this is quantitatively true, we show 
in FigJT] the evolution of the domain length L z along the field direction, for quenches at 
A = 40 and A = 8 in which an external field (a B = 50) was switched on at various times 
during the formation of the bijel. For the larger A, at least, it seems clear that early switch- 
on does in fact promote significantly stronger anisotropy, presumably because the magnetic 
interactions (which result in chain formation) can in part template an anisotropic evolution 
of the interface before jamming is complete. 
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Figure 6: Morphologies for ip = 0.3 with A = 40 at ksT = 2.133 x 10~ 5 ; system size A 3 = 64 3 . 
(a) The initial stable configuration at t = 5 x 10 5 in zero field. The final configurations after 
switching on a field and then running for a further 5 x 10 5 timesteps, with (b) a = 50, (c) 
a = 100, and (d) a = 1000. In panels (e)-(g) are shown the time evolution of characteristic 
length scales of fluid domains in x, y, and z directions; the black dotted line is the initial 
run (a), the red solid line is for case (b), green dashed line (c) and blue dash-dot line is (d). 

14 




Figure 7: The time evolution of the characteristic length scale along the direction of a 
uniform field of strength ag = 50, for various switch-on times tj of the field. Panel (a) is for 
A = 40 and (b) is for A = 8. The black curve is the growth rate of magnetic bijels without 
field; red, green blue, and orange are for t{ = 10 4 , 10 5 and 5 x 10 5 respectively. 

4.2 Field gradient effects 

Our study of the effects of a field gradient is motivated by the experimental work referred to 
already |15j . in which a droplet emulsion stabilized by magnetic particles was controllably 
broken down by exposure to a field gradient. For small field gradients a slight elongation 
of the droplets was seen, but above a critical level the emulsion was completely destabilized 
leading to coalescence of droplets. We plan to address elsewhere [?] the problem of mod- 
elling quantitatively this phenomenon for droplets, but here present data for the analagous 
situation involving bijels. 

First, we generate the initial stable magnetic bijels in a closed-wall simulation box (A 3 = 
64 3 ) quenched at i/jq = 0.3. As usual we set <fi = 0.20. The magnetic colloids again have 
a high dipolar coupling constant A = 40; we use our default temperature setting for which 
the dimensionless parameter describing the relative strength of dipolar and capillary forces 
is A/e = 0.0325. To create a field gradient we introduce an external point magnetic dipole 
exterior to the box, which interacts with the colloidal dipoles; formally this involves addition 
of an extra particle in Eqj3] with a much larger m value than any of the others. This is the 
simplest way to introduce a field gradient and automatically ensures V.B = as required. 
So long as the additional dipole is placed far away from the box, the resulting field gradient 
is (almost) uniform. 

The ratio of field-gradient to capillary forces for a single particle is given by a dimension- 
less number which we now define as C\ = m|V 'B\j-noa. The dipole moments of the external 
magnet are M = 1.64 x 10 6 (Ci = 0.1) and 1.64 x 10 5 [Ci = 0.01) in our simulations, which 
correspond to 5.14 x 10 _11 Am 2 and 5.14 x 10~ 12 Am 2 in the laboratory. The position of the 
external magnet is at (32, 32, -200) for a box size A 3 = 64 3 . We present in FigJH] snapshots 
showing the result of switching on a field gradient for values C\ = 0.01 and 0.1. For the lower 
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value the field gradient only weakly perturbs the bijel, which remains intact. For C\ = 0.1, 
however, the magnetic particles are pulled off the fluid-fluid interface and accumulate at the 
bottom of the container in the region of maximum magnetic field. The interface in the up- 
per part of the container is deprived of particles, and coarsens until further growth becomes 
finite-size limited via interaction with the container walls. 

This destabilization appears closely analagous to that reported by Melle et al [15] for 
droplets. To check whether it actually is similar, we need to compare the relevant dimen- 
sionless control parameter. We first assume this is the single particle force ratio C\ as defined 
above. If as before we identify our simulations with A = 40 as describing (hypothetical) mon- 
odomain magnetite particles of 17nm radius, then the required C\ = 0.1 corresponds in the 
laboratory to a field gradient of roughly 4 x 10 7 T/m. This is not feasible experimentally. 
In fact it is about seven orders of magnitude larger than the field gradient used in the 
experiments of [15]; however, these also used much larger particles (500 nm radius) of a 
superparamagnetic material. Assuming full saturation gives a value C\ ~ 10~ 5 at the point 
where the emulsion droplets are found to break in those experiments. 

Elsewhere we will however show that C\ is not the relevant parameter to describe field- 
gradient-induced particle detachment in complex geometries [?]. This is because the field 
gradient forces, which act like a body force on each particle, generally have a component 
tangential to the interface which is resisted by the interparticle repulsions through the contact 
network. This means that such forces can be transmitted over large distances and can in 
some circumstances act together in expelling a particle from the interface. Indeed, for the 
bijel in a fixed-wall container (FigJBJ) the net force on all particles present is resisted solely 
by the contacts between the bottom layer of particles and the basal wall. Accordingly, each 
of these particles feel a detachment force that is enhanced by a factor N e ff, the number 
of particles that it effectively supports against the field-gradient force, so that the relevant 
control parameter is not in fact C\ but C = C±N e ff. In an undistorted bijel geometry, we 
have N e ff = Nt/Nb — A/2a = 15 where Nt is the total number of particles in the simulation 
and Nb is the number in contact with the base. Although the field-gradient induced breakup 
of a bijel requires unphysical field gradients in the simulations reported here, N e ff increases 
linearly with sample thickness. We cannot simulate very large samples, but if we did, we 
would expect the field gradient to fall accordingly towards the experimentally accessible 
range. Indeed, a laboratory bijel sample of 1 cm thickness, made with particles of 500nm 
radius, might have N e ff ~ 10 4 so that if breakdown occurs at C = 0.1 the threshold value 
of C\ is about 10~ 5 , close to that seen in the experiments of [T5j on droplets. 

5 Discussion 

In this paper, we presented simulations of a binary fluid with magnetic colloids, undergoing a 
quench causing phase separation of the fluids. The neutral wettability of the colloids causes 
them to sequester at the fluid-fluid interface and then jam to form a bijel [HI [TT]. For both 
a symmetric quench and an asymmetric quench, we found the basic morphologies remain 
broadly the same as for the systems with non-magnetic colloids. Dipole-dipole interactions 
cause head-to-tail alignment of the magnetic moments in the plane of the interface, whereas 
uniform external fields align all moments with the field direction. But, unless these interac- 
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Figure 8: Field gradient effects for A = 40 in a closed simulation box. From left to right, 
the initial snapshot at t — 5 x 10 5 ; the configuration after a further 4 x 10 5 timesteps at 
dimensionless field gradient C\ = 0.01 and C\ = 0.1. The magnetic field is vertical. 

tions are simultaneously large, there is little change to either the interfacial particle packing 
or the fluid domain morphology that stems from it. These results show that, although 
dipole-dipole and dipole-field interactions are not negligible in determining the properties of 
magnetic bijels, for typical parameter values the dominant balance is that between capillary 
and short-range repulsive interactions, just as it is for the non-magnetic case. 

Such typical parameter values (such as A = 4 for the ratio of magnetic interactions to 
thermal energies) are those pertaining to magnetic nanoparticles that are already widely 
available for ferrofluid applications. In that context, extremely large A (~ 40, say) may offer 
little advantage since such dipoles will aggregate in a quasi-irreversible fashion. In bijels, 
the magnetic interactions compete primarily with capillary forces, which are much larger 
than thermal ones, so the use of relatively large ferromagnetic particles may be appropriate. 
Moreover, since the structural changes reported here depend on the simultaneous presence 
of large dipolar interactions and large magnetic fields, in practice it may be possible to use 
superparamagnetic particles rather than the permanent dipoles studied in our simulations. 
If so, there is no longer any requirement to have a magnetic monodomain in each particle 
since the paramagnetic particles will anyway achieve saturating magnetization, creating in 
effect larger A than any achievable with monodomain (hence nanometric) particles. Indeed, 
the spectacular experiments of Melle et al. [15] on droplet emulsions use superparamag- 
netic particles of micron dimensions. Note that although it would be possible to perform 
simulations with A > 40, there would be little chance for even local thermal equilibration 
of magnetic degrees of freedom [IB] and so we have not attempted this here. (Indeed for 
nanoparticles the effective run times of our simulations lie in the microsecond range [6j: this 
is long enough for the bijel structure to undergo arrest but not to confirm explicitly that the 
structure is permanent - for confirmation of which we appeal to experiment pTj.) 

When large fields and dipole-dipole interactions were both present, we observed significant 
anisotropy in the fluid domain morphology whose quantitative strength, but not qualitative 
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features, depended on whether the external field was imposed from the start of the quench 
or switched on later. In this study we did not find a combination of dipolar and field 
parameters that would allow the macroscopic morphology to be switched dramatically, e.g. 
from a bicontinuous to a lamellar morphology. It remains possible that by a more complete 
exploration of parameter space such effects could be found. If so this would certainly be of 
interest for applications, allowing the permeability of the bijel to its two component fluids 
to be selectively controlled, in two different directions, by using a magnetic field. 

At modest dipolar interaction strengths we found strong frustration in the head-to-tail 
alignment and a slow relaxation of magnetic degrees of freedom. Even at the longest times 
simulated, significant frustration remained. Based on a careful examination of a single 
droplet arising within an asymmetric quench, we attributed this to an incompatibility be- 
tween global head-to-tail alignment and defects in the particle packing (relative to a hexag- 
onal 2D array) created by the curvature of the interfacial particle layer. These magnetically 
unfavorable defects can be locked in by the strong capillary forces which the dipolar inter- 
actions cannot overcome. 

At first sight, use of a magnetic field gradient (which can exert forces on particles) is 
less promising as a means of morphological control than that of a uniform field (which 
exerts only torques). Indeed the parameter C\ = m\VB\/iTcra, which controls the relative 
importance of field-gradient and capillary forces on a single particle, cannot be bigger than 
the parameter mB/naa 2 which governs the relative importance of magnetic and capillary 
torques in a uniform field. (To achieve equality, the field would have to fall from peak values 
to zero on the scale of one particle diameter.) However, unlike torques, the field-gradient 
forces can be transmitted by contact repulsions within the interfacial layer and can thereby 
gain a cumulative effect. Although in our simulations the resulting amplification factor is 
modest (~ 15), for bijels we can see no limit to it in principle as the system thickness, in 
the direction of the field gradient, is increased. As discussed elsewhere [?] related effects 
probably underly the experimentally field-gradient control of emulsion coalescence reported 
in [15]. It would certainly be interesting to conduct similar experiments on more complex 
interfacial assemblies stabilized by magnetic particles, and we hope that our predictions for 
these and other properties of magnetic bijels might stimulate such experiments in future. 
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